Exact Diagonalisation of The 
XY-Hamiltonian of Open Linear Chains 
with Periodic Couphng Constants and Its 
^^ . Apphcation to Dynamics of 

§: One-Dimensional Spin Systems 

> 

O ■ K. E. Feldman* 

(N 

> ■ DPMMS, CMS, University of Cambridge, Wilberforce Road, Cambridge, CBS OWB, UK 

o ■ 
l> 
o 

' A new method of diagonalisation of the XY-Hamiltonian of inhomogeneous open 

• linear chains with periodic (in space) changing Larmor frequencies and coupling 

[ constants is developed. As an example of application, multiple quantum dynamics 

of an inhomogeneous chain consisting of 1001 spins is investigated. Intensities of 
multiple quantum coherences are calculated for arbitrary inhomogeneous chains in 
^ , the approximation of the next nearest interactions. 

I Key words: linear spin chain, nearest-neighbour approximation, three-diagonal 



X 



Abstract 



matrices, diagonalisation, fermions, multiple-quantum NMR, multiple-quantum co- 
herence, intensities of multiple-quantum coherences. 
PACS numbers: 05.30.-d; 76.20.-hq 

1 Introduction 

One-dimensional exactly solvable models (spin chains, rings) PP have been actively em- 
ployed for studying various problems of spin dynamics [2] and quantum information the- 
ory . Substantial progress in our understanding of spin dynamics has been achieved on 
the basis of a homogeneous XY model for spin chains (s = 1/2) in transverse magnetic 
field |4j. Recently, Hamiltonians of the simplest inhomogeneous systems (alternating sys- 
tems) have been diagonalised for the ring [HI and the linear spin chain [Hj . Development of 
methods for exact solution of inhomogeneous spin problems has become especially urgent 
in conjunction with recent progress in quantum computation and quantum information 
theories j3|. In particular, these methods can be used for studying the quantum state 
transfer from one end of the chain to another one [7| |HI The qubit addressing problem 
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can be attacked by variation of the Larmor frequencies of different spins in classi- 
cal one- dimensional models. This immediately brings inhomogeneity into the XY model 
as diagonal elements of the corresponding Hamiltonian are not equal. Consequently, we 
arrive at one-dimensional spin models with Hamiltonians described by three-diagonal ma- 
trices which elements on the diagonal are not-equal, and those under and over the main 
diagonal are not identically equal as well. Low sensitivity of nuclear magnetic resonance 
(NMR), which is widely applied in experimental implementation of quantum computa- 
tions leads to a further complication of the model described above. Consideration 
of /c-qubit systems brings us naturally to the study of fcn + r spin chains {r < k, n is 
arbitrary, k, n, r are positive integers), with Larmor frequencies and constants of spin-spin 
interactions repeating periodically with period k. 

The paper suggests a new method of diagonalisation of a Hamiltonian of a linear k- 
periodic spin chain of length kn + r for various values r < k. Special attention is given 
to the case k = 3. For a spin chain of 3n + 2 sites with periodic parameters (period 
3) a multiple quantum dynamics 12] is analysed and intensities of all MQ coherences 
are calculated (see Figure 2). This analysis is based on the explicit diagonalisation of 
the Hamiltonian of a linear 3-periodic spin chain with 3n + 2 sites (Theorem 15. 1|) and 
the exact formulae, obtained in the paper (Theorem 16. for the intensities of all MQ 
coherences developed in any nuclear spin system coupled by the nearest neighbour dipolar 
interactions. We conclude the paper with a discussion of the properties of eigenvalues and 
eigenvectors of the general fc-periodic systems with kn — 1 sites f Theorem 17.11) . 



2 Model 

The Hamiltonian of a spin-1/2 open chain with only nearest neighbour (NN) couplings 
has the following general form 

N N-1 

H ^ ^ ^n^nz ~l~ ^ ^ i,In,x-fn+l,x ~l~ In,yln+l,y) ) (1) 

n=l n=l 

where a;„, n = 1, ...,N, are the Larmor frequencies, and Dn^n+i, n = 1, N — 1, are the 
NN coupling constants. The nuclear spins are specified by the spin-1/2 operators Ina at the 
sites n = 1, N with the projections a = x,y, z. The Jordan-Wigner transformation |1] 

-^n,— '^n,x ^^n.y ( 2) ( ]^;_]^ Il^z j C^, 

In,+ = In,x + lln,y = (-2)"-^ {^XllT' hz) C+, (2) 

from the spin-1/2 operators Ina to the creation (annihilation) operators (c„) of the 
spineless fermions takes the Hamiltonian (Q) into the Hamiltonian 

N ^ N~l 

H = ^UJn (C+C„ - 1/2) + - ^ Dn,n+1 (c+C„+i + C++iC„) , (3) 
n=l n=l 
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or in the matrix notations into 



N 



(4) 



ra=l 



In Eq. dD) we denote the row vector (c^, c^) by c+, the column vector (ci, Ctv)* by 
c (the superscript 't' represents the transpose), and specify the matrices fl and D as 



n 



ui 

uji 











ujn-1 

lun 



D 



Di 

Di 

















Dn-i 
Dn-1 



(5) 



To diagonahse the matrix D + 2Q means to construct such unitary matrix U = {unk}nk=i 
and diagonal matrix A = diag{Xi, ...,Xn} that 



D + 2n = UAU^ 



(6) 



where the superscript '+' represents the conjugate transpose, and columns {uu, ...jUniY 
of U, I = 1, N, form an orthonormal basis of eigenvectors of H = D + 2Q. The new 
fermion operators 7^ and 7^ introduced by the relations 



N 



N 



(7) 



k=l 



k=l 



bring Hamiltonian into the Hamiltonian 



1 ^ 1 ^ 



k=l 



n=l 



with energies Aj,/2 of the free fermion waves. 



The problem of diagonalisation of D + 20, with arbitrary constants Un, n = 1, ...,N, 
and Dk, k = 1, N — 1, for large N [N ~ 10^) and computation of various functions of 
D + 2Q is a time consuming problem usually dealt with the help of super computers. On 
the other hand, there are two known exact solutions (suitable for studying systems with 
large N) for such a problem when the spin-system considered has periodic changing Larmor 
frequencies and coupling constants. Namely, the case of equal sites {ui = ... = a;„ = a, 
Di = ... = D]\f_i = c) has been solved in while the case of period 2 with odd (cji = 
uj^ = ... = UJ2 = UJ4 = ... = ujn-1, Di = = ... = Dn-2, = -D4 = ... = -Dat-i) 
has been solved in [H]. The present paper concerns properties of a general periodic chain. 
In the following sections we shall demonstrate how to reduce the diagonalisation problem 
of the Hamiltonian of a fc-periodic chain with kn — 1 sites {k > 2,n > 1) to the problem 
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of finding roots for explicitly given polynomials of degree less or equal k. Therefore, the 
methods developed are particularly useful when n » k. 

In what follows we will utilise essentially the result on diagonalisation of the homoge- 
neous (in other words, 1-periodic) chain jl] which we state now. 

Lemma 2.1 Let c ^ 0, then the matrix 



J„(a,c) 



a 


c 


■■ 


■ 








c 


a 


c ■ ■ 


■ 











c 


a ■ ■ 


■ 














■■ 


a 


c 











■■ 


c 


a 


c 








■■ 


■ 


c 


a 



has n-distinct eigenvalues 



Xj = a + 2c ■ cos 



n + 1 



and the corresponding eigenvectors are of the form 



sm 



sm 



2vrj 

n + 1 J ' V'^ + l 
Proof. It is a straightforward verification. 



j = l,...,n. 



sin 



rnrj 
n + 1 



(9) 



(10) 



3 Reduction over the Period 

Let us consider an open periodic chain of kn + d sites {k > d > 0, k > 1) with fc-periodic 
non-zero NN coupling constants: 



k+l 



D 



kn+l ) 



Dh = D 



k+d 



D 



kn+di 



Dd+1 — Dk+d+i — ... — 
and /c-periodic Larmor frequencies 



n— l)+a!+l ) 



D, = D 



2k 



D 



kni 



(12) 



^1 — ^k+l 



^kn+1, 



^d — ^k+d 



^kn+d-i 



^d+1 — ^k+d+1 — ■■■ — ^k{n-l)+d+l, ' ' ' ^k — ^2k — ■■■ — ^kn, (13) 

{Di,ujj G M, Di 7^ 0, i,j = l,...,k). We shall demonstrate how to reduce the diagonali- 
sation problem for the Hamiltonian of such a system to the problem of diagonalising of a 
certain kxk block matrix which entries are matrices of dimension mxm, (m+1) x (m+ 1), 
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m X (m + 1) and (m + 1) x m {m — n in Section 5 and m = n — 1 in Section 7). A partic- 
ular case of this reduction will be used in the consecutive sections for studying /c-periodic 
chains of length kn — 1. 

Before we proceed we make an agreement regarding our notations. If a Gothic letter 
is used in the description of any matrix as its element, then this means that the corre- 
sponding place in the matrix is a matrix. The size and the structure of a matrix denoted 
by a Gothic letter should be clear or given in the context. Therefore, the original matrix 
is a block matrix. To underline that the matrix consists of blocks we shall often denote 
such matrices in bold. We shall also denote by Im the m x m identity matrix. 

The Hamiltonian of a /c-periodic system with kn + d sites has the following form: 



where 



and 



D 







yj 


n 


n 


n 




n 


n 

yj 


n 


n 







002 ■ 
















• 













■ 















■ 










n 


n 


n 


UJk 


n 

U 




n 

U 


f) 

u 


n 


n 







• 















■ 













• 















• 













• 
















UJl ■ 













■ 
















■ 




1 







• 
















• 










Di 


•• 















•• 








Di 





D2 ■■ 















•• 











D2 


•• 















•• 














•• 





Dk- 


1 







•• 














•• 









Dk 




•• 














•• 





Dk 









•• 














•• 















•• 














•• 















•• 





Dd-1 








•• 















•• 


Dd-1 






(14) 



(15) 



(16) 



To diagonalise matrix H one finds (real) eigenvectors u^, and (real) eigenvalues Ajy, p 
... kn + d, which satisfy the following equation 



{D + 2n)u^ 



= 1, 
17) 
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To resolve this equation we associate to each vector u G M'^"-+'^ vectors 'U(j), j = 1, k, 
formed by those coordinates of u whose numbers have residue j modulo k. Observe that 
among vectors ...,M(a;) there are (i-vectors of dimension (n + 1) and {k — d) vectors 
of dimension n. Now equation ()17|) can be rewritten as a system of n linear equations in 



H 



with matrix H = 2f2 + D. Here 



















= K 




U(k) 









■■■ 





m 


■■■ 






Ws 













(1^ 








>Vfc_i 

m 



(19) 



Wj = ujjIn+1 for j 
degenerate cases 



1, d, yVj = ojjln for j = (i + 1, /c, and with an exception of two 










Vi 



































V 


2 



































V2 





































■■ 



























D = 





■■ 





-Dd 












































Vd+2 









































Vk- 


2 





































Vk-i 








. 'Dk 

























Vk- 


1 




with Vj = 




1 for j 


= 1,2,. 


.,d - 


-1, 


-^^ 




for j 


= rf + 1,.. 


., A; — 1, and 


R"+i - 


R" 


given by 































■■ 


■ 





" 








" 


Dk 





... 












■■ 


■ 




















Dk 


... 













Dd ■ 


. 










Vk = 













■•• 







































■■ 


■ Dd 























... 


D 










■■ 


■ 




. 








. 








... 






(20) 








(21) 
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It is easy to see that 



Dlln, 



VaVi = DaDkJl V^V'^ = D^D^J^, (22) 



where J„ is the Jordan n x n cell: 



1 






1 














1 

1 




(23) 



There are two degenerate cases. In the case k — 2 matrix D decomposes as 



D 



£* 
£ 



(24) 



where JC : 





D2 


•• 


• 













■■ 


• 














Di ■■ 


. 














■■ 














•• 


• 







(25) 



while in the case k > 2 and o? = we do not get matrix Vd and matrix 
of the form 

Dk ••• 
Dk ■■■ 

0' 



IS 










Dk 




(26) 



It is easy to see that this matrix satisfies T>kT>l = Dlln-i,i for a diagonal matrix /n-i,i 
with the first n — 1 diagonal elements equal to 1 and the last one equals to 0. 



Remark 3.1 From the reduction obtained above we divide consideration of the diagonal- 

isation problem into three different cases: period k = 2,3 and k > 4. In each case there 
are further reductions depending on the value of d. In the next sections we shall work out 
the case k > 3 and d = k — 1. 
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4 Some Auxiliary Results 

The diagonalisation process of Sections 5 and 7 will relay on some elementary facts about 



matrices of the form Hij 



H., 



1,3 





j+i where for i < j: 






Di • • • 















Di+i 2a;j+2 • • • 




























D 



3-1 



(27) 



These facts are mostly know and we give their proofs only for the reader's convenience. 

Lemma 4.1 If Dg ^ for all s — i, — 1, then all the eigenvalues of the matrix H^ j 
are distinct. 

Proof. Under the assumption that all D^, s = i, j, are non-zero, one can reconstruct 
every eigenvector Uy of Hi j by its eigenvalue and the first coordinate Ui. Moreover, 
the expressions for all other coordinates -u^,, s = 2, j — z + 1, of arc linear in Ui. This 
implies that every two eigenvectors Uy and u'^ with the same eigenvalue \y are proportional 
with coefficient ui/u[. 

Lemma 4.2 Non-diagonal elements of the adjoint matrix of {Hij — XI 



3-1 



+ij are 



adj,,t {{Hi J - Xlj-i+i)} = adjj^, {{Hij - XIj_i+i)} 



= i-iy^'det (iJ,,,+,_2 - A/,_i) A+.-i ■ ■ ■ D,+t_2det {H,+tj - A/,-.-t+i) (28) 

where the first index in adj^^jj-} denotes the row number, the second denotes the column 
number and t > s. The diagonal elements of the adjoint matrix are 

adj^^^ {{Hij - A/j_i+i)} = det (if^.j+^-i - XIs^i) det (i^i+s+ij - AJj_i__J . (29) 

Proof. The first part follows from an observation that element Di^s-2 situated in row 
(s — 1) and column s of matrix {Hij — A/j_j+i) does not contribute to the cofactor (s, t) 
of {Hi J — A/j_j+i), and similarly for element i^i+t-i situated in row t and column {t + 
1). The second part is due to the splitting of the complement to the element (s, s) of 
{■^i,3 ~ A/j_j+i) into {Hi^i+s-i — Xlg-i) and {Hi+g+ij — Xlj-is)- 

,-1 



CoroUctry 4.1 If we denote elements of {Hij — A/j_j+i) by Ps,t (s — row, t 
umn), then for t > s 



col- 



s.t 



/ ^Y+t ^^t {Hi,i+s-2 — Xlg-i) Dj+s-i • • • Di+t-2det {Hi_^.tj — Xlj^i^t+i) 
^ ' det {Hij - Xlj^i+i) ■ 



The diagonal terms Pt^t df^G 

det {Hi^i+t-i - Xlt-i) det {H,+t+ij 



Pt,t — 



Xlj-i-t) 



det (Hij — Xli 



j-i+i) 



(30) 
(31) 
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Lemma 4.3 If Di ^ 0, then matrices Hi ^ and i?2,fe (k >2) have no common eigenval- 
ues. 



Proof. First, observe that it is true ior k — 2: 

det (ifi,2 - A/2) = (2cui - A) ■ {2uj2 - A) - Dj, 
and if A is a root of det (i/2,2 — A/i) = 0, then A = 2a;2, which is the root of 

{2uji - A) • {2uj2 -X)-Dl^O 



only if Di = 0. 

For general A;, if A is a root of 



det {H2,k - A4_i) = 0, 



(32) 



(33) 



(34) 



then for such A: 



det (i/i,fe - A4) = (2a;i - X)det {H^^k - Xh-i) - Dldet {H^^k - A4-2) = 

= -Dldet {Hs,k - XIk-2) . (35) 

For the last polynomial we can assume that it is non zero by inductive hypothesis and 
because Di ^ 0. From here the statement follows. 



5 Exact Diagonalisation for a Spin Chain with 3n + 2 
Sites 

To demonstrate how the reduction over the period works, first we study a partial case of 
the main result, namely a chain of period 3 with 3n + 2 sites. The reduction of Section 3 
leads us to the following system of linear equations: 



(36) 









" 2Wi 










V 


t 

3 
























2W2 


~ Xi^In+l 


V 


t 
2 






«(2) 


= 0. 














Vs 




V2 




2Wz- 


- Xulj 






. ^(3) . 






Remind that W. , 


3 


= 1,2,3, 


are just scalar matrices, as 


well as the matrix V 






1 




in 


are given 


below: 




















" D2 








... 





" 






' 




•• 


• 










D2 





... 


















Ds ■■ 


• 





V2 = 








L>2 


■-. 


















■■ 


. 










































... D, 


















•• 


• 


Ds 













... 


D2 


. 






. 





•• 


• 






(37) 
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From the second equation of (jHHj) 

(Aj. - 2LJ2) ■ U(2) = ViU^i) + V\u(;i). (38) 

Assume that an eigenvalue Xi, of H is not equal to 2^2, then 

Di , 1 



u 



(2) 



■^(1)+ •^2^(3). (39) 



Substituting ^(2) into the first equation and into the third equation of (jHEl), and using 
that 'D2D\ = D^In we obtain a system: 



(Dl - (A, - 2a;i)(A, - 2^2)) ^^(1) + {DiVi + (A, - 2^2)^^^) ^^(3) = 0, 
((A, - 2uj2)Vs + D{D2) + {Dl - (A, - 20^2) (A, - 20^3)) U(3) = 0. 



(40) 



Lemma 5.1 // 

(A, - 2oj^) (A, - 20J2) ^ Dl (41) 
then the following inequality holds 

(A, - 20^2) (A. - 2073) ^ Dl (42) 

Proof. Observe, that if = 2uj2, then 

(A, - 2c^i) (A, - 2CU2) = Oj^Dl (43) 

and 

(A, - 2CU2) (A, -2iu^) = 0^Dl (44) 
Therefore, it is left to consider the case when Aj, 7^ 2ct;2- Assume that 

(A, - 2002) (A. - 2003) = Dl (45) 

Then from the second equation of 

((A, - 2uj2)V3 + D^V2) = 0. (46) 

With respect to the standard scalar product (-, m = n,n + 1, we have: 

= (m(3), ((A, - 2w2)V3 + DiP2)m(i))r" = {{DiVl + (A, - 2uj2)Vl) M(3), (47) 

Therefore, in order to satisfy the first equation of ()40|1 we necessarily have 

r {Dj - (A, - 2^i)(A, - 2^2)) = 0, 

\ (DiP* + (A,-2cU2)I?*)M(3) =0. ^ 

From here and under the assumption that 

{K-2uj,){X,-2uj2)^ Dl (49) 
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we obtain that = 0. Because Dj 7^ 0, j = 1, 2, 3, the matrix 

{DiVl + (A, - 2lj2)VI) (50) 

has rank n and its kernel is 0. Therefore, ^(3) = 0. Finally, from (jHUj) we obtain that 
U(2) = and thus, Xi, is not an eigenvalue for H. This verifies the lemma. 

Assume as in Lemma (5.11 that 

{X,-2uj,){X,-2uj2)^Dl (51) 
then from the first equation of pUj) we deduce that 

£>iP^ + (A^ - 2iU2) Vj 
{X,-2uj^){X,-2uj2)-Dl 

Substituting into we obtain 



(A.-2^i)P*+AP^ 
= (A.-2..0(A.-2...)-Dr^^^- ^''^ 

Due to Lemma [5. II 

(A, -2^2) (A. -2^3) ^^2- (54) 
Thus, from the second equation of (j40|) we have that 

= (A„-2..)(A„-2a,3)-Bf'- 

Therefore, because X^, 7^ 2ti;2, M(3) satisfies 

((A, - 2cji)(A, - 2cj2)(A, - 2CU3) - (A, - 2uo-i)Dl - (A, - 2uji)DI - (A, - 2uj2)DI) M(3) = 

= {D,VsVl + D,V2Vl)u^,y (56) 
This leads us to the following theorem: 

Theorem 5.1 For j = l,...,n each of the three solutions of the cubic equation 

(A, - 2cui)(A, - 2cu2)(A. - 2CU3) - (A, - 2uj2)Dl - (A, - 2uji)Dl - (A, - 2c^3)D2 = 



= 2D,D2D,cos(^] (57) 

is an eigenvalue for H . The component u (^3) of the corresponding eigenvector is 

f ■ f \ ■ ( '^'^j \ ■ ( ^TTj A A 
= sin , sm , sm . (58) 
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The remaining components and U{2) are determined uniquely fi 



rom 



and 



where T>2, T)^ are given in d^/Tp . 

There are two other eigenvalues of H which satisfy the following quadratic equation 

{K - 2cui) (A, - 2uj2) -Dl = 0. (61) 

The component ^(3) of the corresponding eigenvector Ui, is necessarily zero, the component 
U(i) spans the (one- dimensional) kernel of 

{K - 2uj2)V^ + D1V2 (62) 

and the component U(2) is 

«(2) = 3^-^ ■«(!). (63) 

All eigenvalues constructed are distinct and they exhaust all 3?7. + 2 distinct eigenvalues 
for H . 

Proof. Let us consider the case when 2w2 is an eigenvalue for H and Uy is the corre- 
sponding eigenvector. Then for the components and M(3) we have from the second 
equation of (jSHI): 

W(i) = -i--I?*W(3). (64) 

From the first equation of (IHUj): 

m = [ ■ ^2 + «(3). (65) 

Substituting this into the third equation of (jHE)) we obtain: 

(-^^ ■ V,V\ + -i- { ^^^-'^^^ . V2V\ + V2V'^ + (2C.3 - 20^2)) u 



(3) 



(2^3 - 2^2) - (2^2 - 2^1)^ - ^^^{Jn + 4) ) W(3) = 0. (66) 



Thus, from Lemma f2. II 



2(cu3 - uj2)Dl + 2(cui - uj2)Dl - 2D1D2D3COS { ] = (67) 
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for some j = l,...,n. From here it follows that 2uj2 also satisfies Eq. ()57|) . moreover, 
the formulae for the component M(3) and (as a consequence) for and M(2) coincide 
with those given in the body of the theorem. We obtain that the case A^, = 2uj2 can be 
considered simultaneously with all the other solutions of (jFTj) . 
Therefore, every eigenvalue for H is either a solution for 

(A, - 2uji) (A, - 2002) = Dl (68) 

or it is a solution of ()57p for some j = 1, ...,n. Because H has exactly 3n + 2 eigenvalues 
and all of them are distinct (see Lemma (4 .11) . all the solutions of the quadratic equation 
and of n cubic equations are multiplicity free and pairwise distinct. 

Now the first part of Theorem 15.11 follows from the remarks above, from and 
Lemma Em Indeed if (|^ admits a nontrivial solution for some X^, then either X^, = 2lj2 
or 

/i = ((A, - 2cui)(A, - 2cj2)(A, - 2UJ3) - (A, - 2uj3)Dl - (A, - 2iUi)Dl - (A, - 2uj2)DI) 

(69) 

is an eigenvalue for DiV^V^ + DiV2T>\. The first case has been already discussed, while 
for the second we apply Lemma f2. 11 It follows then that /i is necessarily of the form 

2D,D2DsCos (^^] (70) 

V''^ + 1/ 

for some j = 1, ...,n. As it was explained before, among those 3n eigenvalues constructed 
there are no repeated. The component M(3) of the corresponding eigenvectors are uniquely 
determined by the eigenvalues due to Lemma 12.11 The other two components and 
ti(2) are also uniquely determined by and respectively. This gives us a unique 
eigenvector of H for each solution of (jFTj) . 

The second part of the statement follows from ()40|) because 

{D,Vl + (A, - 2^2)^^3) (71) 

has rank n and its kernel is 0. 



6 Multiple Quantum Spin Dynamics of an Inhomo- 
geneous Spin Chain with 3n + 2 Sites 

Information on the exact spectrum of the Hamiltonian of an open spin chain provides 
us with the techniques for determining the multi-quantum dynamics in such a system. 
The MQ NMR dynamics of the nuclear spins coupled by the nearest neighbour dipolar 
interactions was developed in 0. The corresponding Hamiltonian is 

1 

HmQ = 2 Dn,n+l{In,+In+l,+ + In-In+1-}- (72) 

n=l 
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The Hamiltonian Hmq ()72|) takes exactly the form of the Hamiltonian H (Q) with the Lar- 
mor frequencies u;„ = for all the sites, by making use of the unitary transformation jT2] 

Y = exp {-i-Kh^x) exp {-iirh^x) ■ ■ ■ exp {-mlN-i^x) (73) 

acting on the even sites so that YHmqY^ = H (Eq. ([T)); {un = 0}). The transformation 
Y brings the initial density matrix at the high temperature approximation (see ^T] for 
details): 

N N 

P(0) = E^^-.^ = E (4^^ -1/2) (74) 

to the form 



N 



p(0) = YhY+ = ^(-l)"-i/„_. 



(75) 



n=l 



where we introduce the total polarisation = '^n=i^n,z- The Liouville-von Neumann 
equation {h = 1) 

^^ = [Hmq,p] (76) 

with the Hamiltonian Hmq (f7^ and the initial density matrix p(0) ()74|1 gives us the 
intensities Gn{t) of n = and n = ±2 orders only with the conservation conditions |12[ [TH] 



(77) 



Theorem 6.1 The intensity Gn{t), n = 0, ±2, of MQ coherences of the Hamiltonian 
Hmq ^ are 



Tr [cos' {Hmq ■ t)] Tr [^Au^Hmq ■ t)] 
Go[t) = 77 , G±2[t) = 



Proof. According to [0] 

G2(t)=G_,(t)=i y: e 

with 



2N 



fc=l,3,... n.=2,4,... 



N 



•jk 



lUkie 2 



Ait 



(75 



(79) 



(80) 



where A; are the eigenvalues of Hmq, I = 1, A^, and the unitary matrix U = {uki}ki=i 
diagonalises Hmq- Let us rewrite ()79|) in the matrix form: 

G2it) = G_2(t) = ^Tr (BoAB.A*) , (81) 

where Bq is the diagonal matrix with ones in odd rows and zeroes in even rows, Bi is the 
diagonal matrix with ones in even rows and zeroes in odd rows: 



Bq = diag{l, 0, 1, 0, ...}, Bi = diag{0, 1, 0, 1, ...}; 
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and the matrix A is 



Observe that 
Therefore, 

Using 
and 



A = S{B^-Bo)S*, S = exp{--HMQ-t]. 



{Bi — Bq) Hmq = ~Hmq (Bi — Bq) . 

{Bi -Bo)-S ■ {Bi - Bo) = S* and Tr [{3^] = Tr [{S* 
AA* = A^ = In, a (Si - Bo) A* = {B, - Bo) (S*)' 



Tr [A (5i - Bo) A*] = Tr [{B^ - Bo) A* A] = Tr [{B^ - Bo) AA*] 



(82) 
(83) 

(84) 

(85) 
(86) 



we deduce 



Tr [BoAB.A*] = -Tr [{I - {B, -Bo))A{I + {B^ - Bo)) A*] = 

-Tr [AA* - {B, - Bo) AA* + A{B^- Bo) A* - {B, - Bo) A {B, - Bo) A* 

= -Tr [In - {S*)'] = -\Tr [S' + {S*)' - 21^] = 
4 o 



-Tr 
2 



From here the result follows. 



2i 



^ \ 2 



-Tr [sin^ {Hmq ■ t)] . 



(87) 



Remark 6.1 In ^6] Formula was proposed for alternating spin chains. Theorem \6.1\ 

establishes the same formula for arbitrary spin chains coupled by the nearest neighbour 
dipolar interactions. 



Theorem 16 . II together with Theorem 15.11 allow us to calculate MQ coherence intensities 
of the zero and the second orders, Go{t) and G2{t) + G_2(^), without performing matrix 
multiplication, by solving 0{N) cubic equations (compare with 9J). Let us consider a 
linear spin chain of length 1001 which consists of four-spin fragments represented in Figure 
1. 



D 



12 



D 



23 



D 



34 



Figure 1. The four-spin fragment of the hnear chain. The distances between 
neighbouring spins in the fragments are 2.7, 3 and 3.3 A. The dipolar couphng 
constants are D12 = 27r6096, 1323 = 2ttAAAA, D34 = 27r3339s-i. 
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The distances between neighbouring spins in the fragments are 2.7A, 3A and 3.3A. 
The dipolar couphng constants D12 = 27r6096, D23 = 27r4444, D34 = 27r3339s~^ are used 
in aU numerical calculations. The intensities of MQ coherences for the inhomogeneous 
linear chain with N = 1001 spins consisting of fragments of Figure 1 are shown in Figure 
2. The dynamics behaviour coincides with the one given in p. 

7 The Generalisation of the Method for Spin Chains 
with kn — 1 Sites 



Now we consider the case of a /c-periodic system with kn — 1-sites. According to Section 
3, to diagonalise Hamiltonian of such a system we have to solve the following system of 
linear equations 



(88) 



where 



and 





«(i) 




«(i) 


H 




= K 













H 



i,fe-i 







_ Vk 


• 




2Wi 









2W2 


^2 





V2 


2W3 




















j = 


l,...,k 


-2, 



Vk-i 2Wk 















2Wfc_2 Vk-2 

Vk-2 2Wk-i 



(89) 



(90) 



that 



T^l,k-l = Hi^k-l ® In- 



(91) 



Assume is not an eigenvalue for Hi^k-i (equivalently, is not an eigenvalue for Hi^k-i] 
and consider the matrix 



G 






In-1 



(92) 
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Under the above assumption on A^, the eigenvector equation 











H 




= Xu 






U(k) _ 




_ M(fc) _ 



(93) 



is equivalent to 



G (H — X^Ink-l] 



U(k) 



(94) 



or in another form to 



1 


■■ 


■ 











«(i) 





1 ■■ 


■ 












M(2) 








G 













■■ 


■ 1 












U{k-2) 





■■ 


• 


1 








U(k-1) 




■■ 


• 


Vk-1 


_ 2Wk - \Jn-l _ 






U{k) 



0. 



(95) 



Let us denote elements of the matrix {Hi^k-i — ^uh-i) ^ by Pij {i — row, j — column). 
Using 

{'Hl,k-l — KIn{k-l)) = {Hi^k-l ® In — \uh-l ® In) ^ = (-f^l,fc-l — ^vh-l) ® In, 

we derive from (jH!^ for k — 1 > j > 1 that 

«0) = - {P,,i ■ + P,,k-i ■ Vi_,) (96) 

Substituting this expression with j = 1 and j = k — 1 into the last equation of we 
obtain: 

(Pi,i ■ VJ)i + Pi,fc_iPfcP*_i + Pk-i,i ■ Vk-iVi + Pk-i,k-i ■ I^k-iVi_,) = 

= {2uJk - K)u(k)- (97) 



Because 



Vk-lDl_^ — D^_^In-l, VkV\ — D^In-l, Vk-lD\ — Dk-lDkJn-l, 



VkVi^, = Dk-iDkJn-i, (98) 

with the help of Corollarv l4.1l we derive that if A^, is not an eigenvalue of Hi ^-i, then the 
component U(^k) of the corresponding eigenvector Uu belongs to the kernel of 

M = {det (i/2,fc-i - KIk-2) Dl + det {H^^k~2 - KIk-2) Dl_^- 
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-{2uj, - K)det iH,,k-i - Kh-2)) ■ In-i + i-lfD, ■■■Dk {Jn-i + JLi) ■ (99) 

Since 

{2uJk - Xu)det {Hi^k-i - Kh-2) - det - Kh-2) Dl_^ = det {Hi^k - Kh) 

we can simplify the expression for A4: 

M = {det {H2,k-i - A,4_2) Dl - det {H^^k - Kh) Dl_^) ■ + 

+ {-lfDi---Dk{Jn-l + Ji-l). (100) 

Therefore, 

[det {Hi^k - Kh) Dl-i - det (i?2,fc-i - KIk-2) Dl) "(fc) = 

= i-lfDi ■■■Dk ( J„-i + J*_i) u(^k). (101) 
We are ready to state the main resuh of the paper. 

Theorem 7.1 Each eigenvalue of the Hamiltonian H of a k-periodic system with kn — 1 
sites is either an eigenvalue of Hi^k-i or it is a solution of the equation 

det {H^^k - \uh) - det {H^^k-i - Kh-2) Dl = (-l)^2Di ■ ■ • Dkcos , (102) 

for some j = 1, ...,n — 1. Equation MU^) does not have repeated roots and all k{n — 1) 
solutions constructed from ilU^) are pairwise distinct and are not eigenvalues for Hi^k-i- 
If K is the solution of MU^) for some j = 1, ...,n — 1, then it is an eigenvalue for H 
and the component U(^k) of the corresponding eigenvector Ui, is 

u^k)= (^sin(^^ ,...,sin(^^^-^^y^ . (103) 
The other components U(^j), j = 1, k — 1, are determined uniquely from 

= TT — V ■ [^1 ■ ■ ■ Dj^idet {Hj+i,k-i - Kh-j-i) ■ ^1+ 

+ {-lf-^det (i7i,,_i - A,/,_i) D^--- Dk-2Vi_,] u^k). (104) 

where 'Dk-i, T>k are given in \21]) (d = k ~ 1). Every eigenvalue K of Hi k-i is an 
eigenvalue of H . The component U(^k) of the corresponding eigenvector Uy of H is zero. 
The component U{i) spans the one- dimensional kernel of 

(-If-^D^ ■ ■ ■ Dk-2Vk-i - det {H2,k-i - Kh-2) Vk- (105) 

The remaining components U(^j), j = 2, k — 1, are 

( T^i-iDi ■ ■ ■ Dj^idet (Hj+i k-i — Kh~j~i) /inc\ 
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Proof. The first part of the theorem follows from Lemma 12.11 because H has exactly 
kn — 1 distinct eigenvalues. Indeed, if an eigenvalue X^, of H is not an eigenvalue of 
Hi k_i, then as it was shown above, the component U(^k) of the corresponding eigenvector 

of H satisfies fllOU) . Therefore, satisfies ()102|) for some j = 1, ...,n — 1. 

For the second part of the theorem we again use Lemma I^Tl The component U(^k) of the 
corresponding eigenvector u^, satisfies fllOip . and, therefore, is uniquely determined by X^ 
as stated in Lemma ITTl The remaining components of the eigenvectors u^^) corresponding 
to the solutions of ()l()2j) can be reconstructed from ()9fj|l using Corollarv 14.11 

Finally, for the last part we observe that from Lemma 14.21 and the property of the 
adjoint matrix 



^ det (i/2,fc-i — XyIk-2) 

+ {2uj, - A,) (-1)^' 



det (if2,fc-i — XyIk-2) 
jDi - ■ ■ Djdet {Hj^2,k-i — X„Ik-j-2) _ 



det (iJ2,fc-i — XyIk-2) 

^ det (if2,fe-i — Xylk~2) ^ ^ det (i/2,fc-i — Xi,Ik-2) 

^ ^ adji^j+i {{Hi^k^i - XJk-i)} _ det {Hi^k-i - XJt-i) _ ^ ^^^^-^ 
^ det (i?2,A:-i — XyIk-2) ''' det (-ff2,fc-i — X^Ik-2) 

where k — 1 > j > 1 and 6ij is the Kronecker symbol. We also have 

Dk-z{ i^k-3 Di---Dk-3i'2uJk-i- X^) _^ (^2ujk-i - X^) ( Di---Dk-2 ^ Q 

det (i/2,fc-i — XyIk-2) det (i/2,fc-i — XyIk-2) 

(108) 

If X^ is an eigenvalue for Hi ^^i, then 

g^.^ I jj^^ y Pidet {Hs^k-i - XJk-2) ^ det{Hi^k-i - XJk-2) ^ ^ ^^^g^ 
det (iJ2,fc-i — XuIk-2) det {H2,k-i — XyIk-2) 

Therefore, if spans the kernel of 

(-If-^D^ ■ ■ ■ Dk-2Vk-i - det {H2,k^, - XJk-2) (110) 



and U{k) = 0, then vector u^, with the other components given by ()106p does satisfy the 
eigenvalue equation 

iH-Xhn.^)u, = 0. (111) 

This completes the proof. 
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8 Conclusion 

In the paper we proposed a new one- dimensional exactly solvable model for a linear k- 
periodic (in space) open spin chain with kn — 1 sites. For the diagonalisation procedure 
it is important that the number of sites is {k — l)(mod/c). Nevertheless, the model can 
serve as a good approximation for any linear periodic spin system if the number of sites 
in it is much more than the period. 

The developed method of diagonalisation of the XY-Hamiltonian of inhomogeneous 
linear spin chains can be applied to different problems of quantum information theory 
[3,7,8] and spin dynamics. This method allows us to avoid matrix multiplications which 
are time consuming operations in the systems with large numbers of spins. In some cases 
we can suggest analytical methods for problems of spin dynamics instead of the known 
numerical ones [9]. 

The suggested method could also be applied to different physical and technical prob- 
lems which use three-diagonal matrices ^3]. In particular, new numerical methods of 
solving such problems could be worked out on the basis of the approach proposed in the 
paper. 
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